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Abstract 

We derive optimal upper and lower bounds on the slope of the Isgur-Wise 
function at the zero recoil point in terms of the sum of the TBB couplings, 
estimated recently from a QCD Spectral Sum Rule (QSSR). The problem 
is solved by means of a duality theorem in functional optimization. Optimal 
correlations between the slope and the convexity parameters of the Isgur-Wise 
function with the same input are also obtained. 



1 Introduction 



In the last time, there has been some interest in the investigation of the 6-number 
form factor F(q 2 ) of the 5-meson 

< B(p')\V»\B(p) >= {p + pTF(q 2 ) (1.1) 

where V 1 = Iry^b and q = p' — p. In the large quark mass limit, this form factor 
coincides with the renormalized Isgur- Wise function [1] £(a>), uj = v-v' = 1— g 2 /2M|, 
describing the semileptonic decays of B into D and D* mesons. The values of the 
form factor in the physical region relevant for these decays {uj > 1, i.e. q 2 < 
) are of interest for the determination of the mixing matrix element \V c b\ from 
experimental data. Various models for the Isgur- Wise function have been proposed 
in the last time in the frame of the heavy-quark effective theory (HQET). On the 
other hand, the possibility of obtaining model independent bounds on the values of 
F(q 2 ), particularly on its slope at the origin (the charge radius) was outlined for the 
first time in [2] . These bounds follow from analyticity, unitarity and the asymptotic 
perturbative QCD expansion for the amplitude of the vacuum polarization due to 
the current. As noticed in the subsequent papers [3-6], the derivation [2] suffered 
from the drawback that it did not include the three upsilon poles present in the 
form factor F below the BB threshold. The problem was correctly treated in [7,8], 
where bounds on F(q 2 ) were derived using as input only the position of the three 
upsilon poles below the BB threshold. Although the methods are different, the 
results obtained in [7] and [8] are both optimal. With a weak form of the unitarity 
[7], and using as input only the poles location, known from the experimental masses 
of the T resonances, the following bounds on the slope of the Isgur- Wise function 
were obtained 

-6. < F'(l) < 4.1, (1.2) 

where the derivative is with respect to uj. The upper bound is actually not interest- 
ing, since it is much larger than the Bjorken upper bound -0.25 [2]. In [7] bounds on 
F'(l) were also calculated by assuming that the residua of T the poles are known. 
Using as input several particular values for the TBB couplings, the conservative 
lower bound 

F'(l) > -1.5 (1.3) 

was proposed in [7], but the authors stress upon the fact that the only rigorous lower 
bound on the slope of the 6-number form factor at the origin is (1.2). 

It is of interest to know whether this bound can be improved using only model 
independent information about the TBB couplings. Recently, using a QCD Spectral 
Sum Rule (QSSR), Narison [9] obtained an estimate for the sum of the residua 
of the three upsilon resonances T(IS'), T(2S'), T(3S*) situated below the threshold 
4M§. However, in [9] this estimate was used only as an indication in choosing several 
particular values for the TBB couplings. Upper and lower bounds on -F'(l) were 
then computed with the method presented in [7], for these specific values of the 
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couplings. No attempt was made up to now to evaluate the real constraining power 
of the relation derived in [9], i.e. to calculate bounds on the slope of the Isgur-Wise 
function in terms of the sum of the upsilon pole residua. 

In the present paper we address this specific problem. We derive upper and lower 
bounds on the slope F'(l) using as input the positions of the T poles and the sum 
of their residua. We also derive an optimal inequality relating the slope and the 
convexity parameters of the Isgur-Wise function, using the same input. These are 
the best model independent theoretical results which can obtained without choosing 
particular values for the TBB couplings. 

The paper is organized as follows: in the next Section we derive upper and 
lower bounds on the slope ^'(1), using analyticity, unitarity, the QCD one-loop 
expansion for the polarization function and the sum of the TBB couplings derived 
from QSSR. We solve the problem by applying a duality theorem in functional 
optimization [10,11]. We also generalize the method to include higher derivatives of 
the form factor at the origin. The results, expressed as explicit inequalities involving 
the quantities of interest, are investigated numerically in Section 3. 



2 Derivation of the bound 

We consider the 6-number form factor F(q 2 ) defined in (1.1), which is a real-analytic 
function in the complex plane t = q 2 , cut along the real axis from 4M^ to infin- 
ity, except for iV = 3 poles below the BB threshold, corresponding to the states 
T(IS'), T(2S*) and T(3S). Their contribution to the form factor can be parametrized 
as [7] 

F(t) = F(0) + t± ^* BBh < + (2.1) 

i=l m Ti 1 lt 

where jt i denote the coupling constants which govern the electronic widths of the 
Tj-resonances 

4-7T 

r(T^ (7) - e + e~) = f 2 Mr -a 2 (2.2) 

and g-x % BB are the coupling constants of the Tj resonances to the BB system. As in 
[7] we denote by rji the product of the above coupling constants 

Vi = ZgriB&hv i = 1, 2, 3. (2.3) 

By using a QSSR for the vertex function, Narison [9] estimated recently the sum of 
these couplings: 

3 

5>* = a = -0.224. (2.4) 

i=i 

We shall take this model independent relation as input in deriving bounds on the 
derivative of the Isgur-Wise function at the origin. We use also the normalization 
condition 

F(0) = 1 (2.5) 
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and the integral inequality [7] 



' r {V 1)3/2 \F(4M 2 B y)\ 2 dy < 1 (2.6) 



12 16Mlx(Q 2 ) h y3/ 2 (y + Q 2 /AMl) 2 
where y = ^r- Here 

m = jVc. f 1 x 2 (l-x) 2 dx 

W j 2W0 m 2 + Q 2 x(l-x) 10 

is the QCD one-loop level expression of the derivative of the polarization function, 
for spacelike values Q 2 > 0, and N c is the number of colours. In deriving (2.6) it was 
assumed [7] that the unitarity sum for the polarization function is saturated with the 
lowest BB states and that the contribution of the intermediate states B + B~, B°B° 
and B®B® is the same in the limit where the light quark mass differences are small. 
The stronger version of unitarity used in [2,6,8], which requires additional dynamical 
assumptions, beyond HQET, will not be adopted below. 

In the inequality (2.6) the spacelike momentum Q 2 enters as a parameter. The 
dependence of the results on the choice of Q 2 , as well as on the QCD sum rule 
(dispersion relation) written for the polarization function was discussed in [8]. In 
what follows we shall take the particular value Q 2 = 0. Then (2.6) becomes 



5n f 



/oo 
(y _ l)3/2 y -7/2 \ F{y)l 2 dy < L ^ 



16iV c 

We write this inequality in a canonical form by performing the conformal mapping 



z = ^IH- 1 (2.9) 

which applies the cut y-plane onto the interior of the unit disk \z\ < 1, such that 
z(0) = and the beginning of the cut, y — 1 is transformed into z = — 1. In the 
new variable the positions of the T poles are denoted by Zi and the residua of the 
form factor F(z) at these poles can be written as 

1 — Zi 

n = (z- Zi)F(z)\ z=Zi = riiZi——^ (2.10) 

1 + Zi 

with rji defined in (2.3). Also, the inequality (2.8) becomes 

-jf w{9)\F{9)\ 2 d6 <1, (2.11) 

where 

W (Q) = ^l cos ^ sin - (2.12) 

v ; 48 2 2 y 1 
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is a nonnegative weight. Of interest for us are the slope and the convexity parameters 
of the form factor 

P 2 = -F'(l) , c = \f"{1) . (2.13) 

the derivatives being with respect to u>, which are related to the derivatives of F 
with respect to z through 

F' z (0) = -8p 2 , F' z '(0) = 128c - 32p 2 . (2.14) 

We shall obtain upper and lower bounds on the derivatives of F, using as input 
the analyticity of the form factor in the cut t-plane, except for the T poles, and the 
conditions (2.4), (2.5) and (2.6) ( the last one being written alternatively as (2.11)). 

In order to bring the problem to a standard form, we consider first the function 
[2,7] 

4>{z) =0(O)(l + ^) 2 v / r^, 0(0) = ^^^, (2.15) 

which is analytic and nonzero inside the unit disk \z\ < 1 and whose modulus on 
the boundary is related to the weight w(9) given in (2.12) through 

|0(0)| = (2.16) 

Let us define now the new function 

g(z) = F(z)<j>(z)B(z) (2.17) 

where F is the form factor and 



D[z) = ^hrk- (2 - 18) 

i=i 1 Zi z 



is a product of so-called Blaschke factors [10] (in our case, to a good approximation, 
z* = Zi). It is clear that the function g is analytic in \z\ < 1, the poles of F being 
canceled by the corresponding zeros of B. Moreover, the constraints (2.4), (2.5) and 
(2.11) on the form factor can be written immediately in terms of the function g. 
Thus, using (2.10) and (2.15) we have 



g{zi) = ri4>{zi 
and the condition (2.4) becomes 



B(z) 



Z - Zi 



(2.19) 



£7tf(*i)=a- (2-20) 



i=l 
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Here the weights 7, are defined as 



7» 



with 



1 

0(0)^(1 + * i )(l-* i ) 3 / 2 £ i 



2 - -2,; 



As concerns the conditions (2.5) and (2.11), they become 

,7(0) = 0(0)5(0) 



and, respectively 



MU 2 



f 1 

I2W0 ' 5 



! d0 



1/2 



< 1, 



(2.21) 



(2.22) 



(2.23) 



(2.24) 



where we took into account (2.16) and the fact that the Blaschke factors have modu- 
lus equal to 1 along the unit circle. Using a standard terminology [10], (2.24) means 
that the function g belongs to the unit sphere of the Hilbert space H 2 . 

We express also the derivatives of g(z) at the origin in terms of the derivatives 
of F(z), which are further related to the slope and convexity parameters through 
(2.14). Using the relations 



0'(O) = ^0(0), 0"(O) = 



(2.25) 



which follow from (2.15), we obtain 



g'(0) = 0(0)5(0) 



3 
2 



B'(0) 
5(0) 



g"(0) = 0(0)5(0) F"(0) + F'(0) 



3 + 2 



B'(0) 
5(0) 



1 , 0^(0) 
" 4 + 6 5(0) 



+ 



5"(0) 
5(0) 



(2.26) 



where we took into account (2.5) and the derivatives of 5 follow from (2.18). 

We pass now to the calculation of bounds on the derivatives of the analytic 
function g in the origin, using as input the constraints (2.20), (2.23) and (2.24). We 
first consider only the slope of the form factor, i.e. the derivative g'(0). As will be 
clear, the generalization to higher derivatives is straightforward. 

In order to solve the problem we shall apply a method based on a norm mini- 
mization over a convex set in H 2 and a duality theorem in functional optimization 
[10,11]. A similar technique was applied in [12] for other problems in particle physics. 
More precisely, let us consider the following set of functions: 



K = \geH 2 \ g(0), g'(0) given, J>^) = a. 



(2.27) 
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It is clear that K is a convex set in the Hilbert space H 2 . In defining it, we took into 
account the constraint (2.20) and the fact that g(0) is known, according to (2.23). 
In addition, we ascribed to the derivative g'(0), which is the quantity of interest, a 
definite value, for the moment arbitrary. 

We have still to impose the condition (2.24) and this will restrict the allowed 
range of the unknown parameter g'((0). More precisely, if this parameter is correctly 
chosen, all the functions in the set K will have L 2 norms less than or equal to 1. 
If, on the other hand, g'(0) is outside the allowed range, the L 2 norms of all the 
functions in K will be greater than 1, and this will be true also for the smallest of 
these norms. It follows that the inequality 



where the l.h.s. denotes the minimal L 2 norm over the set K, describes the optimal 
range of the parameter g'(0). As we shall see, the minimum norm (2.28) will depend 
explicitely on the specific information used in the definition of K, which consists 
from the values g(0),g'(0) and the real number a. This will allow us to obtain 
explicitely rigorous upper and lower bounds on the unknown value g'{0) in terms of 
the value g(0) and of the parameter a. 

Having reduced the problem to the calculation of a minimum norm, we resort 
for solving it to a duality theorem in functional analysis [10,11]. Usually, a duality 
theorem relates a minimization over an abstract space to a maximization in the dual 
space, which is often simpler than the original one. More precisely, if K is a convex 
set in the Hilbert space H 2 and g denotes an arbitrary function belonging to K, the 
theorem of interest for us has the form [10,11]: 



Here the supremum in the right hand side is calculated upon all the complex func- 
tions G defined on the boundary of the unit circle, of L 2 norms less than 1. In 
writing (2.29) we took into account the fact that the dual of the space L 2 is again 



We have first to calculate the supremum over g e K in the r.h.s. of this equa- 
tion. In a standard terminology [10] this supremum is called the support functional 
of the convex K. Of course, of interest are only the finite values of the support func- 
tional. This clearly restricts the form of the the functions G(z) taken into account 
in the functional optimization. As we mentioned above, the function G entering the 
maximization (2.29) is a complex function of L 2 norm less than 1, being in general 
the boundary value of a function having an analytic part inside the unit disk, and 
possible interior singularities. By residua theorem, these singularities pick up in the 
integral (2.29) the values of g at some interior points. It is clear that these values 
must be known, according to the definition of the set K, otherwise the supremum 



min \\g\\ L 2 < 1, 



(2.28) 




(2.29) 



L 2 . 
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upon g G K would be uncontrolably large. By applying this argument to the partic- 
ular set defined above we see that the functions G yielding a finite support functional 
of K have the general form 



G(z) 



01 _|_ ° 2 _|_ ^2 



z 



1=1 



Z - Zi 



+ Q(z) 



(2.30) 



where the residua ai,a 2 , bi and the analytic function Q(z) are arbitrary. By inserting 
this expression in (2.29) we obtain 



mm \\g\\L 2 = sup sup 

geK GeH 2 ,\\G\\ L 2<igtK 



a ig (0) + a 2 g'{0) + J2big{zi) 

■i=i 



(2.31) 



where we applied the residua theorem, which implies in particular that the analytic 
function Q does not contribute. 

It is clear that the supremum with respect to g in the r.h.s. of this equation is 
finite only if the residua bi are proportional to the weigths 7$ entering the constraint 
in the definition of K. We must have therefore 

bi = a 3li (2.32) 

where a 3 is an arbitrary parameter, along with ai and a 2 . Then (2.31) becomes 

min \\g\\ L 2 = sup[aig>(0) + a 2 g'(0) + a 3 a], (2.33) 

where the parameters entering the maximization in the r.h.s. must satisfy the 
condition that the L 2 norm of G is less than 1. By performing the Fourier analysis 
of the function G of the form (2.30) we obtain after a straightforward calculation 
the quadratic constraint 



2 Hlkzfzl 



(ai + a 3 7i) 2 + (a* + a 3 J2 7*^) 2 + + a i T 

i=l i=l i,k=l 



< 1, 



(2.34) 



which must be satisfied by the parameters a*. 

The linear optimization problem (2.33) with the quadratic constraint (2.34) can 
be solved immediately, by means of a linear change of variable which brings this 
constraint to a diagonal form. The solution can be written in a compact form, and 
using (2.28) we finally obtain the inequality 



^(0)+^(0) + - 



-, 2' 



a-s(o)X>-</(o)£7* 



where 



A 



i=i 



3 22 



i=l 



< 1, 



El! f« 



:,fc=i 



ZiZk 



(2.35) 



(2.36) 
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Using (2.23) and (2.26) the inequality (2.35) can be written in terms of F(0) and 
F'(0), allowing one to calculate upper and lower bounds on the slope of the Isgur- 
Wise function using as input the sum , a, of the T pole residua. The numerical 
results are presented in Section 3. 

Before closing this section we mention that the technique applied above can be 
generalized in a straightforward way as to include higher derivatives of the form 
factor in the origin. This can be done by a suitable definition of the convex set K 
and a proper choice of the function G used in the duality theorem (2.29). Instead 
of (2.27) we now define 



K = \geH 2 \ g(0), g'(0) , g"(0) given, = «• 



i=i 



and the function G of (2.30) is replaced by 



(2.37) 



G(z) = — + — + -3 

z z 2 Z 6 



E 

1=1 



Z - Zi 



+ Q(z). 



(2.38) 



A calculation similar to the one presented above leads to the following inequality 
involving the first two derivatives of the form factor: 



where 



«-s(o)E7,-s'(o)E7, 



e=Ef 

i,k=l 



i=l 



-1 2 s 



Y~ E iiA 



i=i 



< 1, 

(2.39) 
(2.40) 



The inequality (2.39), written in terms of the slope and convexity of the Isgur-Wise 
function by means of (2. 23), (2. 26) and (2.14), yields a rigorous correlation between 
these parameters and the sum a, of the T poles residua. We shall apply it in the 
next Section for testing specific models of the HQET. 



3 Numerical results and conclusions 

In our analysis we took M B = 5.3 GeV and varied m h in the range 4.7 — 5.3 GeV. 
We first calculate upper and lower bounds on the derivative -F'(l), using (2.35). 
The results are presented in Table 1, for several values of the parameter a in the 
constraint (2.4). 
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a 


F'(l) ■ 


F'(l) 


-2.240 


-3.05 


4.30 


-0.272 


-3.91 


3.77 


-0.224 


-4.13 


3.61 


-0.0224 


-4.23 


3.55 


0. 


-4.23 


3.53 


0.224 


-4.36 


3.45 


0.672 


-4.55 


3.29 


2.240 


-5.15 


2.59 



Table 1 : Upper and lower bounds on the slope of the Isgur-Wise 
function for various values of the sum (2.4) of the T pole 
residua. 



In Table 1 we investigated for completeness several values of a, up to ten times 
larger than the value proposed in [9] and given in (2.4). Values with opposite sign 
were considered as well. Taking a = —0.224 we obtain 

-4.13 < F'(l) < 3.61 (3.1) 

which slightly improve the values (1.2) and is the best bound on the slope of the 
6-number form factor, obtained using only model independent information about 
the residues. 

We consider now the inequality (2.39) relating the slope and the convexity pa- 
rameters of the Isgur-Wise function to the sum of the TBB couplings. According 
to (2.39), the domain of the allowed values for these parameters is the interior of 
an ellipsa in the plane p 2 — c. It is of interest to establish the compatibility of 
some particular models of HQET with this rigorous inequality. For illustration, we 
consider the following parametrizations proposed in [9]: 

F(u) ~ 1 + F'(l)(w-1), 

F(u) ~ exp[F'(i)(u-l)}, 

r i 1 2F '« 

F( W )~{-(l+a;)} , 

F(u>) ~ -^-exp |[2F'(1) + . (3.2) 

i -+- ou i uj -+- 1 ) 

These parametrizations, with the choice F'(l) ~ —1.0, were used in [9] for a new de- 
termination of the mixing angle V c b- They correspond, respectively, to the following 
expressions for the convexity parameters: 

c = 0. 
c=\[F\l)f 
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c=-F'(l)[2F'(l)-l] 



c — - 



F\lf 



F'iD - - A 



(3.3) 



In order to test the compatibility of the models (3.2) with the constraint (2.4), we 
calculate the derivatives g'(0) and g"(0) using the relations (2.14) and (2.26) and 
we insert them in (2.39). The result of this calculation is given in Table 2, which 
indicates the l.h.s. of (2.39) denoted by LHS, for several v alues of p 2 and the 
parameter c computed according to (3.3) for each specific model (3.2). The two 
columns correspond to a — —0.224 as proposed in [9], and also to a value three 
times larger, respectively. Unlike the bounds on the slope given in Table 1, which 
had a slow variation with a, the results are now quite sensitive to the sum of the 
pole residua. According to (2.39), a result less than 1 is expected to occur, if the 
parametrizations (3.2) are consistent with the constraint (2.4). However, as seen 
from Table 2, for p 2 = 1 a value greater than 1 is obtained for all the models (3.2), 
which shows that these models are not compatible with the relation (2.4) proved 
in [9]. This throws some doubt on the determination of V c b made in [9], which was 
based on the models (3.2) with the choice p 2 = 1. On the other hand, for p 2 = 0.5, 
the models (3.2) turn out to be consistent with the relation (2.4), while they fail to 
be so if the value of a is increased by a factor of 3. 



P 2 


c 


LHS 


LHS 






a = -0.224 


a = -0.672 


1. 


0. 


1.65 


1.95 


1. 


0.5 


1.37 


1.78 


1. 


0.75 


1.39 


1.83 


1. 


0.875 


1.44 


1.89 


0.5 


0. 


0.91 


1.29 


0.5 


0.125 


0.87 


1.28 


0.5 


0.25 


0.86 


1.30 


0.5 


0.25 


0.86 


1.30 



Table 2 : Left hand side of (2.39), denoted as LHC, for p 2 = 1 
and p 2 = 0.5, and c computed from (3.3) for the models 
(3.2). The two columns correspond to a = —0.224 and 
to a value three times larger. 



In conclusion, in the present paper we succeeded to fully exploit the constraining 
power of the QSSR relation (2.4) between the TBB couplings, derived in [9]. The 
problem was solved by means of some powerful techniques in functional optimization. 

The results show that the relation (2.4) has a rather weak constraining power 
on the slope of the Isgur-Wise function. In particular, it slightly improved the lower 
bound on F'(l) from the value -6, when no information on the TBB couplings 



10 



is used, to the value -4.13, which is still very generous. On the other hand, the 
constraint (2.4) is able to produce nontrivial correlations between the slope and the 
convexity parameters of the form factor. The relation (2.39) derived in the present 
paper turns out to be an useful tool for investigating the compatibility of the specific 
models of the quark theory with analyticity, unitarity and QCD. 
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